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Abstract 

We study in this paper nonlinear subdivision schemes in a multivariate set- 
ting allowing arbitrary dilation matrix. We investigate the convergence of 
such iterative process to some limit function. Our analysis is based on some 
conditions on the contractivity of the associated scheme for the differences. 
In particular, we show the regularity of the limit function, in and Sobolev 
spaces. 

Keywords: Nonlinear subdivision scheme, convergence of subdivision 
schemes, box splines 
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1. Introduction 



Subdivision schemes have been the subject of active research in recent 
years. In such algorithms, discrete data are recursively generated from coarse 
to fine by means of local rules. When the local rules are independent of the 
data, the underlying refinement process is linear. This case is extensively 
studied in literature. The convergence of this process and the existence of the 

n n 

limit function was studied in 2!] and [7i] when the scales are dyadic. When the 
scales are related to a dilation matrix M, the convergence to a limit function 

n n n 

in was studied in [8| and generalized to Sobolev spaces in [5] and [13j. In 
the linear case, the stability is a consequence of the smoothness of the limit 
function. 

The nonlinearity arises naturally when one needs to adapt locally the re- 
finement rules to the data such as in image or geometry processing. Nonlin- 
ear subdivision schemes based on dyadic scales were originally introduced by 
Harten 9| 10| through the so-called essentially non-oscillatory (ENO) meth- 
ods. These methods have recently been adapted to image processing into 
essentially non-oscillatory edge adapted (ENO-EA) methods. Different ver- 
sions of ENO methods exist either based on polynomial interpolation as in 
Qjll or in a wavelet framework jsl, corresponding to interpolatory or non- 
interpolatory subdivision schemes respectively. 

In the present paper, we study nonlinear subdivision schemes associated 
to dilation matrix M. After recalling the definitions on nonlinear subdivi- 
sion schemes in that context, we give sufficient conditions for convergence in 
Sobolev and L'p spaces. 



3 



2. General Setting 

2.1. Notations 

Before we start, let us introduce some notations that will be used through- 
out the paper. We denote the cardinal of the set Q. For a multi-index 

fj, = (/ii, fj,2, ■ ■ ■ , fid) G and a vector x = {xi, X2, - ■ ■ , Xd) G M*^ we define 

d d d 

M = 1^- = n l^i- and xi^ = n 

i=l 1=1 i=l 

For two multi-index m, e N'^ we also define 

mi J y 

Let ^(Z'^) be the space of all sequences indexed by Z'^. The subspace of 
bounded sequences is denoted by i'^{Z'^) and ||M||^oo(2d) is the supremum of 
{\uk\ : k e Z'^}. We denote the subspace of all sequences with finite 

support (i.e. the number of non-zero components of a sequence is finite). 
As usual, let &{'L'^) be the Banach space of sequences u on Z*^ such that 
[[^^[[^^(zd) < where 

llwIl^pCz-i) := X] l^fcl^ for 1 < p < oo. 




As in the discrete case, we denote by L^(M'^) the space of all measurable 
functions / such that ||/||LP(Md) < oo, where 



LPi^d) ■■= (^J ^ \f{x)\Pdx^ for 1 < p < oo 



and ||/||L°°(Md) is the essential supremum of |/| on E". 

Let /i G N'^ be a multi-index, wc define the difference operator 
■ ■ '^d^j where V^^ is the fXjth difference operator with respect to the 
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jth coordinate of the canonical basis. We define as D^^ • • ■ D'^'^, wfiere 
Dj is the differential operator with respect to the jth coordinate of the canon- 
ical basis. Similarly, for a vector a; e M'' the differential operator with respect 
to X is denoted by D^. 

A matrix M is called a dilation matrix if it has integer entries and if 
lim M~" = 0. In the following, the invertible dilation matrix is always 

n—^oo 

denoted by M and m stands for \det{M)\. 

For a dilation matrix M and any arbitrary function $ we put $^^^(2;) = 

We also recall that a compactly supported function $ is called I^-stable 
if there exist two constants Ci, C2 > satisfying 

Ci||c||£p(zd) < II ^ Cfc$(x - A;)||ip(]Rd) < C2||c||^p(zci). 

Finally, for two positive quantities A and B depending on a set of param- 
eters, the relation A < B imphes the existence of a positive constant C, 
independent of the parameters, such that A < CB. Also A B means 
A < B Si^d. B < A. 

2.2. Local, Bounded and Data Dependent Subdivision Operators, Uniform 
Convergence Definition 
In the sequel, wc will consider the general class of local, bounded and 
data dependent subdivision operators which are defined as follows: 

Definition 1. For v e ^^{'L'^), a local, hounded and data dependent subdivi- 
sion operator is defined by 

S{v)wk = ^ak-Mi{v)wi, (1) 
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for any w in £°°(Z'^) and where the real coefficients ak^Mi{v) G K are such 
that 

ak-Mi{v) = 0, if \\k - Ml\\eoc(j^d) > K (2) 

for a fixed constant K . The coefficients ak{v) are assumed to he uniformly 
hounded hy a constant C , i.e. there is C > independent of v such that: 

\ak{v)\<C. 

Note that the definition of the coefficients depends on some sequence v, while 
S{v) acts on the sequence w. 

Note also that, from and (j2]) the new defined value S{v)wk depends 
only on those values I satisfying \\k — M/||£oo(2d) > K. The suhdivision oper- 
ator in this sense is local. 

To simplify, in what follows a data dependent subdivision operator is an 
operator in the sense of Definition [H With this definition, the associated 
subdivision scheme is the recursive action of the data dependent rule Sv = 
S{v)v on an initial set of data v^, according to: 

= Sv^-^ = S{v^-^y-\ J > 1. (3) 

2. 3. Polynomial Reproduction for Data Dependent Suhdivision Operators 

The study of the convergence of data dependent subdivision operators 
will involve the polynomial reproduction property. We recall the definition 
of the space Pat of polynomials of total degree A^: 

P^:={P;P(x)= ('f^^'}- 
With these notations, the polynomial reproduction properties read: 
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Definition 2. Let N > be a fixed integer. 

1. The data dependent subdivision operator S has the property of reproduc- 
tion of polynomials of total degree N if for all u G £°°(Z'^) and P E^n 
there exists P G Pat with P — P E Pat-i such that S{u)p = p where p 
and p are defined by pk = P{k) and pk = P{M~^k). 

2. The data dependent subdivision operator S has the property of exact 
reproduction of polynomials of total degree N if for all u G £°°(Z'^) 
and P G Pat, S{u)p = p where p and p are defined by pk = P{k) and 
Pk = P{M-'k). 

Remark: 

The case = is the so-called "constant reproduction property". For 
a data dependent subdivision operator defined as in ([1]), the constant repro- 
duction property reads Yl (^k-Mi{v) = 1, for all v G £°°(Z'^). 

3. Definition of Schemes for the Differences 

Another ingredient for our study is the schemes for the differences asso- 
ciated to the data dependent subdivision operator. The existence of schemes 
for the differences is obtained by using the polynomial reproduction property 
of the data dependent subdivision operator. 

Let us denote A' = (V^, = I) and then state the following result on 
the existence of schemes for the differences: 

Proposition 1. Let S be a data dependent subdivision operator which re- 
produces polynomials up to total degree N . Then for 1 < I < N + 1 there 
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exists a data dependent subdivision rule Si with the property that for all v,w 

in t^ilJ^), 

A'5(v)w := Si{v)^^w 

Proof: 

Let I be an integer such that 1 < Z < + 1. By using the definition of 
with = I, we write: 

Prom the definition of S{v)w we infer that 



max(//i,---,/id) 



V^S{v)Wk = (-1)' I ^ I XI 0,k-m-e-Mp{v)Wp, 

mi,— ,md=0 \ ^ 



where we have used the notation m - e — m\ei H — • + maed- Straightforward 
computations give 

max(^ti,- - ,^td) f fJ' \ 

V^S{v)Wk = X^Wp X (-l)M \ ak-m-e-Mp{v) 



mi,— ,md=0 \ 
= ^Wpfk,p{v,il). (4) 



Let us clarify the definition of fk,p{v, /i). Since the data dependent subdivision 
operator is local we have ak-Mp{v) = for any data v e £°°[Z'^) and any index 
k such that \\k — Mp||^oo(2d) > K. Now by putting k — e + Mn, we get that 
fk,piv, IJ>) is defined for p in the set 

V^'ik) := {p: ||n-p + M-^(£-m -6)1100 < K\\M-^\\^, < < Hi^i} 

Then, we define V{k) := {p : \\k — Mp||oo < K}. Since the data dependent 
subdivision scheme reproduces polynomials up to total degree N, we have 



for any |z/| = r < N: 

J2 ak-MpW = P,{k) forallfcGZ^ (5) 

p&V{k) 

where is a polynomial of total degree r. By tacking the differences of 
order \u'\ = r + 1 in ([5]) we get 

Note that the above equality is true for any u such that |z/| = r. We deduce 
that {fk,p{v,i''))i^ G Z'^ is orthogonal to {p'')p^v''' (k) "where |g| < r. Note that 
\^(ySn-i3)nev-'{k) ,W\ = r + 1, (3 e spans {p'^)p^v-' (k) ^^^^ 
write for any p G V^' (k): 

\u\=r+l ■r6Z<* 

Now, by using (jl]) we obtain for any \fi\ < N + 1: 

peVt^ik) \u\=lr€Zd 
p£Vi^{k) \i/\=l 

If we now make fi vary, we obtain the desired relation. □ 

Now that we have proved the existence of schemes for the differences, we 
introduce the notion of joint spectral radius for these schemes, which is a 



generalization of the one dimensional case which can be found in 
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Definition 3. Let S{v) : P'iJJ^^ — P'iJJ^) he a data dependent subdivision 
operator such that the difference operators Si{v) : {P'{l/)Y'' — ?■ (£^(Z'^))'", 



with qi = |/i| = /} exists for I < + 1. Then, to each operator Si, 

/ = 0, ■ ■ ■ , + 1 (putting Sq = S) we can associate the joint spectral radius 
given by 

1 

PpAS) := inf ||(5^)^||;p(2d),r 
In other words, Pp{S) is the infimum of all p > such that for all v G 

£P(Z'^), one has 

\\A^S^v\\t,p(^zd),, < p^||A'i;||£p(2;d),,, (6) 

for all j > 0. 

Remark: Let us define a set of vectors {xi, ■ ■ ■ , x^} such that [xi, ■ ■ ■ , a;n]Z" = 
(i.e. a set such that the linear combinations of its elements with 
coefficients in Z spans Z'^). We use the bold notation in the definition of 
the set so as to avoid the confusion with the coordinates of vector x. Then, 
consider the differences in the directions xi, ■ ■ ■ ,Xn- One can show that there 
exists a scheme for that differences which we call Si for / < + 1 provided 
the data dependent subdivision operator reproduces polynomials up to de- 
gree N (the proof is similar to that using the canonical directions). If we 
denote by A' the difference operator of order / in the directions xi, ■ ■ ■ ,Xn, 
one can see that || A^wH^p^^d)?; ~ || A'i;||fp(2d)'3i for all v in £p(Z'^) and where 
Qi = #{/^! I/^I = h = (/^i)i=i,- - ,n}-Then, following ([6]), one can deduce that 
the joint spectral radius of Si is the same as that of 5*^. 

4. Convergence in L''' spaces 

In the following, we study the convergence of data dependent subdivision 
schemes in W which corresponds to the following definition: 
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Definition 4. The subdivision scheme = Sv^~^ converges in Lp(M'^), if 
for every set of initial control points G P'ilJ^), there exists a non-trivial 
function v in L^(R°'), called the limit function, such that 

lim \\vj - vWivmd) = 0. 

d 

where Vj{x) = Yl '"i4>j,k{x) with = Yl niax(0, 1 — \xi\). 

4.I. Convergence in the Linear Case 

When 5* is independent of v, the rule ([1]) defines a linear subdivision 
scheme: 

SVk = ^ ttk-AIlVi. 

If the hnear subdivision scheme converges for any v G P'{'L'^) to some function 
in U'iW'') and if there exists such that hm ^ 0, then {a^, k G Z*^} 
determines a unique continuous compactly supported function $ satisfying 

= ^ afc$(Mx - k) and ^ $(x - A;) = 1. 
Moreover, t>(x) = ^ f°$(a; — A;). 

4-2. Convergence of Nonlinear Subdivision Schemes in U Spaces 

In the sequel, we give a sufficient condition for the convergence of non- 
linear subdivision schemes in Lp(M°'). This result will be a generalization of 
the existing result in the linear context established in and only uses the 
operator 5*1. 

Theorem 1. Let S be a data dependent subdivision operator that reproduces 
the constants. If Pp^i{S) < mv, then Sv^ converges to a L^ limit function. 
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Proof: Let us consider 



d 

where 4>{x) = J| max(0, 1 — \xi\) is the hat function. With this choice, one 

i=l 

can easily check that ^ 4>{x — k) — 1. Let Pp,i{S) < p < m^, it follows that 



||A^'t;^||(^p(zd))ci < p^'||A^'t;°||(^p(zci))d, 
since qi — d. We now show that the sequence vj is a Cauchy sequence in I/: 

Vj+i{x) - Vj{x) = j2 ^rVj+i,fe(2;) - j2 ^p</'j,p(^) 
= Yl zi(^r^ -'f^p)0j+i,fe(^)</'j,p(^) 

where we have used J2 0(" ~ ^) — 1- Now, since the subdivision operator 
reproduces the constants: 

peZ'i A;ez<* lez^ 

Note that 



Since E afc_MK^')-<^p-; = 0, {V,5,_;3, ^ e {F(k) U {p}} , ^ e i = 1, 

spans {ttk-Mi — 5p-i)i&{F{k)vj{p]}- This enables us to write: 

d 

Vj+i{x) - Vj{x) = Y Y^k,p,i^yi<l>j+i,ki^)'PjA^)- 

pez<« feezd i&v{k) UM «=i 
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Since | ^ (f)j^i^k{x)\ = 1 following the same argument as in Theorem 3.2 of 
, we may write: 



\vj+i-Vj\\LP(Rd) < m p max \\ViV^\\ep(zd) 

l<i<d 



~ (Ayi|Ait;°||,.(z.) (8) 

which proves that Vj converges in L^, since p < rrip . Note that, for p = oo, 
we obtain that the limit function is continuous. □ 

Furthermore, the above proof is valid for any function $o satisfying the prop- 
erty of partition of unity when p = oo. In general, we could show, following 
Theorem 3.4 of ^|, that the limit function in is independent of the choice 
of a continuous and compactly supported $o- 

4-3. Uniform Convergence of the Subdivision Schemes to C functions (s < 
1) 

We are now ready to establish a sufficient condition for the smoothness 
of the limit function with s < 1. 

Theorem 2. Let S{v) be a data dependent subdivision operator which re- 
produces the constants. If the scheme for the differences satisfies Pp^i{S) < 
m~^^p , for some < s < 1 then Sv^ is convergent in LP and the limit 
function is . 

i_ 

Proof: First, the convergence in is a consequence of pp^i{S) < nip . 
In order to prove that the limit function v be in C"^, it suffices to evaluate 
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\v{x) —v{y) \ for — y||oo < 1- Let j be such that m ^ ^ < \\x — y\\oo < ""^ "'• 
We then write : 

\v{x)-v{y)\ < \v{x)-Vj{x)\ + \v{y)-Vj{y)\ + \vj{x)-Vj{y)\ 
< 2||t; - + \vj{x) - Vj{y)\ 

Note that (jS]) imphes that \\v — Vj\\x^oc(^d) < p-' || A^t>°||£oo(2d). Since Vj is 
absolutely continuous, it is almost everywhere differentiable, so putting y = 
X + M~^h, with h = (/ij)j=i,...,d satisfying ||/i||oo < 1 we get: 

\vj{x + M~^h) -Vj{x)\ < \vj{x + M'^h) -Vj{x + M~^{h- hded))\ 

+ \vj{x + M-\h - hded)) - Vj{x + M-\h - hdCd - hd^ied-i))\ 
+ ■■■ + \vj{x + M-^{hiei)) - Vj{x)\ 

Then, using a Taylor expansion we remark that, there exists 9 e] ~ hd,hd[ 
such that: 

\vj{x + M-^h)-Vj{x + M-'{h-hded))\ = ^v{hdDd(j){M^x + h-k + eded) 

If we denote ^d{,x) = $(xi) ■ ■ ■ ^{xd-i)'^{xd), where is the characteristic 
function of [0, 1] and ^(xj) = max(0, 1 — we may write: 

\v,{x + M-^h)-Vj{x + M-\h-hded))\ ~ ^ Vdi;^/id^d(y) 

where yi = {M~^ x + h — k + 9ded)i ifi < d and yd = 2{M~^x + h + 9ded)d — kd 
(we have used the fact that the differential of the hat function $ is the Haar 
wavelet). Iterating the procedure for other differences in the sum, we get: 

d 

\vj{x + M~^h) - Vj{x)\ < ^ ||Vif^||£^(zd) < ||A^f^||foo(2d). 
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Combining these results we may finally write: 

\v{x)-v{y)\ < \v{x) - Vj{x)\ + \v{y) - Vj{y)\ + \vj{x) - Vj 

< 2||^; - Vj\\L^(^^d) + \vj{x) - Vj{y)\ 

< piAV||^oo(z.) + ||AV 

< ^■||AV||,oc(z.) < 



loo 



with s < - log(poo,i)/ logm. □ 

5. Examples of Bidimensional Subdivision Schemes 

In the first part of this section, we construct an interpolatory subdivision 
scheme having the dilation matrix the hexagonal matrix which is: 



M 



2 1 
-2 



For the hexagonal dilation matrix, the coset vectors are eq = (0, 0)^, ei = 
(1, 0)^, £2 = (1, -1)^, ss = (2, -1)^. The coset vector ei,i = 0, • • • , 3 of M 
defines a partition of as follows: 

3 

= [j{Mk + ei,k eZ^} . 

The discrete data at the level j, is defined on the grid P = M~^l?, 
the value vl. is then associated to the location M~^k. We now define our bi- 
dimensional interpolatory subdivision scheme based on the data dependent 
subdivision operator which acts from the coarse grid V^~^ to the fine grid 
grid To this end, we will compute at the different coset points on the 
fine grid using the existing values v^~^ of the coarse grid T^~^, as follows: 

15 



for the first coset vector Eq — (0,0)-^ we simply put I'lffe+ep = ^fie 
coset vectors i — 1, • • • ,3. tfie value v;^^,.,.^^, i = 1, • • • , 3 is defined by 
afiinc interpolation of the values on the coarse grid. To do so, we define four 
different stencils on T^~^ as follows: 

Vi'' = {M-^+'k, M-^+\k + e,),M-^+\k + 62)}, 

V^'^ = {M-'+^k, M-'+\k + 62), M-^+\k + 61 + 62)}, 

W^'^ = {M-^+\k + 6i),M-^+\k + 62),M-^+\k + 6i + 62)}, 

Wl = {M-^+^k,M-^+\k + 6i),M-^+\k + 6i + e2)}. 

We determine to which stencils each point of belongs to, and we then 
define the prediction as its barycentric coordinates. Since we use an affine 
interpolant we have: 

1 ■ 1 • 

^iik = and v^Mk+e, = 2^r^ + 2^fe+er (9) 

To compute the rules for the coset point £2, Vk <^an be used leading 

respectively to: 



'"Mk+e2 ~ ~'~ ^^k+e2 ^^k+ei+e2- ^^^) 



When one considers the rules for the coset point £3, Wl or can be used 
leading respectively to: 

"^Mk+es ~ ^^k+e-i + ^^k+ei+cz + 2 

This nonlinear scheme is converging in L°° since we have the following result: 
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Proposition 2. The prediction defined by l[W\) . IpJ]) satisfies: 

i\ A \ j II ^IIaI? '111 

II A 'yM.+eJu°°{Zd) < ^11 A V-^ hoc(^zd) 

We do not detail the proof here but the result is obtained by computing 
the differences in the canonical directions at each coset points. If we then 
use Theorem [2] we can find the regularity of the corresponding limit function 
is C with s < -^^^^^^ ^ 0.207. 

log (4) 

In the second pat of the section, we build an example of bidimensional 
subdivision scheme based on the same philosophy but this time using as 
dilation matrix the quincunx matrix defined by: 



M 



-1 1 
1 1 



whose coset vectors are Eq = (0,0)"^ and ei = (0,1)^. Note that ao,o = 1 
and since the nonlinear subdivision operator reproduces the constants we 
have (^Mi+e = 1 for all coset vectors e. To build the subdivision operator, 

i 

we consider the subdivision rules based on interpolation by of first degree 
polynomials on the grid ^Mfc+ei corresponds to a point inside the cell 

delimited by M~^~^^{k, k + ci, k + 62, k + ci + 62}. There are four potential 
stencils, leading in this case only to two subdivision rules: 



V 



Mk+ei 



Vr + Vile.^e,) (12) 



VMk+e, = iK+ei +^^fc+ej (13) 



Note also, that since the scheme is interpolatory we have the relation: v-' 



Mk 



vl ^. Let us now prove a contraction property for the above scheme. 
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2. when k = Mk' + ei, we can show that: 



I ill ^ II A 1 h 



Proposition 3. The nonlinear subdivision scheme defined by fT^) and / f7g|) 

satisfies the following property: 

1. when k = Mk' : 

II i.l i II ^ II A 1 7-2 II 



1 

2' 

The proof of this theorem is obtained computing all the potential differ- 
ences. This theorem shows that the nonlinear subdivision scheme converges 
in L°° since pi,oo('S') < 1. 

6. Convergence in Sobolev Spaces 

In this section, we extend the result established in [is!] on the convergence 
of linear subdivision scheme to our nonlinear setting. We will first recall the 



notion of convergence in Sobolev spaces in the linear case. Following [12 1 
Theorem 4.2, when ^o{x) = Yl (^k^oi^x — k) is L^'-stable, the so-called 
"moment condition of order k + 1 for a" is equivalent to the polynomial 
reproduction property of polynomial of total degree k for the subdivision 
scheme associated to a. In what follows, we will say that $0 reproduces 
polynomial of total degree k. When the subdivision associated to a exactly 
reproduces polynomials, we will say that $0 exactly reproduces polynomi- 
als. We then have the following definition for the convergence of subdivision 
schemes in Sobolev spaces in the linear case 



13|: 
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Definition 5. We say thatv^ — Sv^ ^ converges in the Sobolev space W^{W^) 
if there exists a function v in W^{R'^) satisfying: 

lim ll^^j - 'y|Ufc(Kd) =0 

where v is in W^{W^), and Vj — ^ vl^o{M^x — k) for any $o reproducing 

feez^ 

polynomials of total degree k. 

We are going to sec that in the nonhnear case, to ensure the convergence 
we are obhged to make a restriction on the choice of $o- We will first give 
some results when the matrix M is an isotropic dilation matrix, we will 
also emphasize a particular class of isotropic matrices, very useful in image 
processing. 

6.1. Definitions and Preliminary Results 

Definition 6. We say that a matrix M is isotropic if it is similar to the 
diagonal matrix diag{ai, . . . , Ud), i.e. there exists an invertible matrix A such 
that 

M = A^^diag{ai, ad) A, 
with \ai \ — . . . — \a(i\ being the eigenvalues of matrix M. 

Evidently, for an isotropic matrix holds \ai\ — ... — \ad\ — cr = m^. 
Moreover, for any given norm in W^, any integer n and any v we have 

(j^lkll < llM^-ull < (7"|k||. 

II II rvj II II rsj II II 

A particular class of isotropic matrices is when there exists a set ei, 62, • • • , 
such that: 

Mci = Aie^(i) (14) 
19 



where 7 is a permutation of {1, • • • ,q}. Such matrices are particular cases of 
isotropic matrices since — XI where / is the identity matrix and where 

d 

X = Yl Xi. For instance, when d = 2, the quincunx (resp. hexagonal) matrix 

•i=i 

satisfies = 21 (resp. = 41). 

We establish the following property on joint spectral radii that will be 
useful when dealing with the convergence in Sobolev spaces. 

Proposition 4. Assume that S reproduces polynomials up to total degree N. 
Then, 

for all n — 0, . . . , N . 

Remark: If M is an isotropic matrix and S reproduces polynomials up to 
total degree A^, then 

Pp,n+l{S) > a~^Pp^n{S), 

for all n = 0, . . . , AT. 

Proof: It is enough to prove 

According to the definition of spectral radius there exists p > Pp,i{S) such 
that for any 

Using the notation cu^ := S{u^^^) ■ . . . ■ S{u^^)u wc obtain 

||Vcu-'||£p(zd) < p^\\Vu\\^p(^zd)- 
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Since 

n 

where 

Kn = (^l-Ml,-r{u^~^)ai^.i-Mlj.2 " • • • " ah-Mn{u°)- 

We can write down the £^-norm as follows: 



where {ej}™^ are the representatives of cosets of First note that: 

||A;-n||oo < \\k~n + M-^ei\\oo + \\M~^e{\\oo. 

Note that M~^el belongs to the unit square so that ||M~%^||oo < Ki. When 
A-' . i ^ 0, one can prove that there exists Ko > such that 

\\k-n + M-^e]\\^ < Ko, 



the proof being similar to that of Lemma 2 in ll|. From these inequalities 



it follows that if A' . ^0 there exists K^i > such that 

\\k - n\\oo < K3, 

that is, for a fixed k, the values of uj for / G {M^k + ej}^^ depend only on 

Un with n : {||A; - n\\oo < K^}. 

Let us now fix k and define u such that 



Ul 



ui, if ||A; - /||oo < -ft'a; 



0, otherwise. 
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Let := S{u^-^) ■ S{vP)u, then 



^0, if - M-^Z||oo > i^4, 
since ii Aj^y^ 0, then 

\\k - M-H\\^ < \\k - n|U + ||n - M-H\\^ < + K2 := K^. 

Moreover, from ||fc-M--'7||oo < K^, it follows that ||M^A;-/||oo < K4M^\ 
Taking all this into account, we get 

E E Ki' = E E i-/i'^E E 

le{Mik+si} k&d i^{Mik+ei} keZ^ \\Mjk-l\\oo<K4Mj 

< ||M||i,||A^a;/||,.(^.) < {\\M\Uy\\A'u 



That is, ||u;-''||^F(zd) < (||M||ooP)-''|K||^F(zd), consequently ^^(5") < ||M||ooP. 
Now, if p ^ pp,i{S) we get pp{S) < \\M\\^pp,i{S). □ 

6.2. Convergence in Soholev Spaces When M is Isotropic 

First, Let us recall that the Sobolev norm on iy^(M'^) is defined by: 



LPiRd) + ^ ||-D^/||LP(R<i). (15) 
\H\<N 



If one considers a set Xi, ■ ■ ■ ,Xn such that [x-^, - ■ ■ , XnjZ" = Z'^, an equivalent 
norm is given by: 



LP{Rd) + ^ \\D''f\\iP(j^dy (16) 



where D'' = Dt^^ ■ ■ ■ Dti- 



We then enounce a convergence theorem for general isotropic matrix M: 
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Theorem 3. Let S be a data dependent nonlinear subdivision scheme which 
exactly reproduces polynomials up to total degree N — 1, then the subdivision 
scheme Sv^ converges in W^{M.'^), provided $o is compactly supported and 
exactly reproduces polynomials up to total degree N — 1 and 

1 s 

Pp,n{S) < mp^'^ for some s > N. (17) 
Proof: Note that because of Proposition HI the hypotheses of Theorem 

1 1 s-l 

[3]imply that Pp,k{S) < nid pp^k+i{S) < rrip , which means that (fl7j) is also 
true for k < N. Let us now show that Vj is a Cauchy sequence in L^. To do 
so, let us define 

d 

1=1 

where A = is defined in For a multi-index /x = (/xi, . . . G Z'^ 
let 

Since A is invertible, the set {g^ : = A^} forms a basis of the space of all 
polynomials of exact degree A^, which proves that 

Now, we use the fact that, since M is isotropic, q^{D){f{M^x)) = a^^{q^{D)f){M^x) 

d n 

where cr'^ = H '^i' We can thus write: 

i=l 

q,{D){v,^, - V,) = q,{D) ( t;/+^<|.o(M^+V - /) -J^^IMM'x - /) ) . 
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We use now the scaling equation of $o to get 



- E <9k-Mr)q,{D) ($o(M^+^x - /)) 
J2 J2iai-Mriv^) - 9i-MrKa''^'^'\q,iD)%)iM^+'x - I). 



Since S and $o exactly reproduce polynomials up to total degree — 1, we 
have for \fi\ < N — 1: 



Remark that gi-Mr = for ||/ — Mr\\ > K since $0 is compactly supported. 
Since ^VSi.p, \u\ = N,r E F{1) = |||/ - Mr\\ < max{K , K)^ , /3 E Z'^j spans 
{ai-Mr{v^) - gi-Mr)r&Fii), we dcducc: 

'. r&F{l) \u\=N 



Consequently, 

Since Pp,n{S) < v(i^I'p~^I'^ ^ with s > we obtain 

\\q^{B){v,^^ - t;,)llLP(M^) < a^'(^-^)l|A^t;°||(,.(^.)),,. 

From this we deduce that —v^\i^v(^d^ tends to with j. Making 

\i vary, we deduce the convergence in iy^(M'^) □ 

We now show that when the matrix M satisfies ( fT4l) and when $0 is a box 
spline satisfying certain properties, the limit function is in W^{^^\ Before 
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that, we need to recall the definition of box splines and some properties that 
we will use. Let us define a set of n vectors, not necessarily distinct: 

Xn = {Xi,--- C Z'^\{0}. 

We assume that d vectors of X„ are linearly independent. Let us rearrange 
the family X„ such that = {x-^, ■ ■ ■ ,Xd} are linearly independent. We 

d 

denote by [xi, ■ ■ ■ , Xd] [0, If'^ the collection of linear combination s V! X iX\ with 



Xi G [0, 1[. Then, we define multivariate box splines as follows |4||16|: 

if X e [xi, - ■■ ,Xd][0,l['^ 



/3o{x,Xd) 



det(xi,--- ,x^)\ 

otherwise 



-1 



Po{x,Xk) = / /3o(x-tXk,Xfc_i)dt, n>k>d. (18) 



'0 



One can check by induction that the support of /3o(x, Xn) is [xi, X2, • • • , Xn ] [0, 1] 



The regularity of box splines is then given by the following theorem 16|] : 



Proposition 5. /3o(x, X„) is r times continuously differentiable if all subsets 
of Xn obtained by deleting r + 1 vectors spans M.'^. 

We recall a property on the directional derivatives of box splines, which 
we use in the convergence theorem that follows: 

Proposition 6. Assume that X„ \ Xj. spans R"', and consider the following 
box spline function s{x) = ^ Cfc/3o(x — fc, X„) then the directional derivative 
of s in the direction Xy. reads: 

DxrSix) = ^ V^.^Cfc/3o(x - fc, X„ \ Xr). 



We will also need the property of polynomial reproduction which is 
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16|: 



Proposition 7. ///3o(x,X„) is r times continuously differentiable then, for 
any polynomial c{x) of total degree d < r + 1, 

p(a;) = J]c(z)/3o(a;-z,X„) (19) 

is a polynomial with total degree d, with the same leading coefficients (i.e. the 
coefficients corresponding to degree d). Conversely, for any polynomial p, it 
satisfies ( fTPI) with c being a polynomial having the same leading coefficients 
as p. 

Theorem 4. Let S he a data dependent nonlinear subdivision scheme which 
reproduces polynomials up to total degree N — 1 and assume that M satisfies 
relation [14\ ), then the subdivision scheme Sv^ converges in W^{M.'^), if when 
N > 2, is a C^^"^ box spline generated by satisfying ^o{x) = 

J2 gk^o{Mx-k) and if N = 1 ^o{x) = ^ gk^o{Mx-k) and ^ ^o{x-k) = 

k k fceZd 

1 and if 

1 s 

Pp,n{S) < mp^'^ for some s > N. (20) 

Proof: We here prove the case N > 2, the case = 1 can be proved 
similarly. First note that since ^o{x) is a C^^"^ box spline, we can write for 
any polynomial p of total degree — 1 at most: 

p{M-^x) = ^p(0<l>o(M-ix-z,Xj 

= XI 9q-MiP{i)^o{x - q, Xn) 

Using Proposition [7] we get p and p have the same leading coefficients, and 
that ^ gq-Mip{i) is a polynomial evaluated in M~^i having the same leading 
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coefficients as p. That is to say the subdivision scheme {Sv^)g = Qq-Mi^l 
reproduces polynomials up to degree N — 1. 

As already noticed, the joint spectral radius of difference operator is in- 
dependent of the choice of the directions Xi,--- , Xn that spans Z*^. Fur- 
thermore, it is shown in IJ], that the existence of a scaling equation for $0 
implies that the vectors Xi, i = 1, ■ ■ ■ , n satisfy a relation of type (IT^ . We 
consider such a set {xi}i=i,...,n and then define $o(a;) = /3o{x,Yn) the box 
spline associated to the set 

{N N 
•^li ' ' ' 1 Xl, ■ ■ ■ , Xyii ■ ■ ■ , X][ 

which is by definition. We then define the differential operator D^^j_j := 
^M-^xi ' ' ' ^M-ix ■ ^^^^ characterization f|T6|) of Sobolev spaces 

therefore = (/ii)i=i,... ,n- For any < we may write: 

-EE vi9p-MiiD^Po)iM^+'x - p, Fat), 



using the scaling property satisfied by /3o- Then, we get: 

where Yj^ is obtained by removing /Xj vector Xi, z = 1, ■ ■ ■ , li to Yat and 
= (V^^!) ^. As both ttk-M.i^'') 9k-M. reproduce polynomials up 
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to total degree N — 1, there exist a finite sequence Ck,p such that: 

^^{J2(^k-Mi{v^) - gk-Mi)vi) = Yl E 

*GZ<i p&v{k)\jvik)W\=M 
where V{k) = {i, \\k - Mi\\ < K}, where gk-Mi = if - Mi\\ > K. We 
finally deduce that: 

D'-^.UvM^) - vj{x)) = E E E c,,,(^^)VX/3o(M^+^a: - k, F^). 
From this, we conclude that: 

Now, consider a sufficiently differentiable function / and remark that DM-j-^xifi^) 
{Df){x).M~^~^Xi, where Df is the differential of the function /. We also 
note that = XI which imphes that X — a'^ and we then put j + 1 — 
q X ['^J + r with r < q and where |_.J denotes the integer part. Prom this 
we may write: 

Dm-^-^.J{x) = (J-^L^J(L'/)(x).M-^xi 

and then 

where rj depends on j. Making the same reasoning for any order /i of differ- 
entiation and any direction Xi, we get, in L^: 

We may thus conclude that 
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To state the above result, we have used the fact that the joint spectral radius 

is independent of the directions used for its computation. Since we have the 

-— - 

hypothesis that Pp^\^\{S) < mp for s > we get that 



which tends to zero with j, and thus the limit function is in W^{M.'^) M. 
A comparison between Theorem 12] and S] shows that when the subdivision 
scheme reproduces exactly polynomials, which is the case of interpolatory 
subdivision schemes, the convergence is ensured provided $o also exactly re- 
produces polynomials. When the subdivision scheme only reproduces poly- 
nomial the convergence is ensured provided that $o is a box spline. Note also 
that the condition on the joint spectral radius is the same. We are currently 
investigating illustrative examples which involve the adaptation of the local 
averaging subdivision scheme proposed in jol] to our non-separable context. 

7. Conclusion 

We have addressed the issue of the definition of nonlinear subdivision 
schemes associated to isotropic dilation matrix M. After the definition of 
the convergence concept of such operators, we have studied the convergence 
of these subdivision schemes in and in Sobolev spaces. Based on the study 
of the joint spectral radius of these operators, we have exhibited sufficient 
conditions for the convergence of the proposed subdivision schemes. This 
study has also brought into light the importance of an appropriate choice of 
$0 to define the limit function. In that context, box splines functions have 
shown to be a very interesting tool. 
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